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INTRODUCTION 

Manual  palpation  remains  the  first  diagnostic  line  of  defense  against  breast 
cancer.  Whether  performed  in  a  self-examination  or  clinical  physical  exam,  breast 
palpation  is  a  low  cost,  totally  non-invasive  test  that  often  is  the  initial  means  of  disease 
detection  [1-4].  Unfortunately,  physical  examination  of  the  breast  has  many  limitations 
including  variable  skill  of  the  examiner  (particularly  for  breast  self-exams)  and 
subjectivity.  Perhaps  its  greatest  limitation,  however,  is  that  sensitivity  to  palpate  masses 
deep  in  the  breast,  or  within  relatively  dense,  stiff,  or  heterogeneous  tissues,  is  poor.  The 
principal  objective  of  this  project  is  to  develop  a  means  to  measure/image  tissue  elasticity 
(i.e.,  “hardness”)  deep  within  tissues.  The  hypothesis  is  that  remote  measurement  of 
elasticity  in  breast  tissues  is  possible  and  provides  unique  information,  which  could 
increase  detection  and/or  characterization  of  potentially  malignant  masses  not  accessible 
to  manual  palpation.  Our  preliminary  studies  suggest  that  the  proposed  methods  are 
capable  of  precisely  measuring  intemd  deformation  and  strain  in  three  dimensions. 
These  data  are  required  to  reconstruct  the  elasticity  distribution  within  the  object. 
Consequently,  technologies  developed  within  the  scope  of  this  project  may  have 
significant  diagnostic  value  in  detection  and  management  of  breast  cancer. 

Imaging  scientists  have  attempted  to  electronically  extend  the  touch  of  the 
physical  examiner  by  a  variety  of  techniques  that  infer  tissue  elasticity.  The  essential 
element  of  these  approaches  is  to  measure  relative  internal  motion  and  strain  of  tissue 
structures  experiencing  mechanical  stress.  Most  elastography  to  date  has  utilized 
ultrasound  [5-14],  although  MRI  has  also  been  used  [15-22].  Usually  an  external  static  or 
dynamic  deformation  is  applied  while  the  resultant  displacement  or  propagating  shear 
wave  is  documented  by  imaging  devices.  In  the  original  MRI  method,  spatial 
magnetization  tagging  was  applied,  but  this  had  poor  spatial  sampling  of  displacement 
and  only  measured  2-D  motion  [18].  More  recently,  motion  phase  encoding  via  bi-polar 
gradients  was  employed  to  produce  2-D  strain  maps  in  media  mechanically  driven  by 
external  forces  [19-22].  This  approach  is  applicable  to  measure  quasi-static  displacement 
or  dynamically  record  propagation  of  mechanical  shear  waves  and  produce  strain  images. 
While  strain  depends  on  internal  elastic  properties,  it  is  also  strongly  affected  by  the 
applied  deformational  geometry,  and  thus  is  not  uniquely  reflective  of  tissue  properties 
(i.e..  Young’s  modulus).  “Dynamic”  measurements  exhibit  local  shear  wavelength, 
which  allow  direct  calculation  of  Young’s  modulus.  Shear  wave  attenuation,  interference 
from  standing  waves  off  multiple  reflectors,  and  limited  resolvable  points  over  the  shear 
wavelength  are  significant  challenges  being  address  by  investigators  of  these  methods. 

Relative  to  ultrasound,  nuclear  magnetic  resonance  (NMR)  imaging  (MRI)  has 
two  major  advantages  for  elasticity  imaging.  First,  the  resolution  and  accuracy  of  NMR 
displacement  and  strain  components  can  be  isotropic  compared  to  ultrasound  estimates. 
Indeed,  ultrasound  can  accurately  measure  axial  (i.e.,  along  the  beam  axis)  motion  at  high 
spatial  resolution  (<mm),  but  lateral  displacement  is  measurable  at  much  lower  spatial 
resolution  defined  by  the  acoustic  beam  width  which  varies  with  depth.  The  third 
dimension  is  generally  not  even  considered  given  the  limitations  of  ultrasound. 
Consequently,  the  “elastogram”  is  compromised  by  reduced  motion  dimensionality  and 
overall  low  motion  resolution  of  the  imaging  system.  In  MRI,  isotropic  resolution  of 
«lmm  can  be  readily  achieved  by  standard  means.  Secondly,  NMR  displacement  and 
strain  estimates  are  less  sensitive  to  the  strain  decorrelation  effects  that  produce 
substantial  error  in  ultrasound  cross-correlation  speckle  motion  estimates.  Indeed,  large 
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strains  change  the  spatial  distribution  of  the  speckle  patterns  used  for  motion  tracking 
markedly  increasing  errors  in  displacement  and  strain  estimates  [13,  14].  This  would 
have  a  significant  impact  on  the  signal-to-noise  (SNR)  of  the  NMR  displacement,  strain, 
and,  subsequently,  elasticity  images. 

Limitations  of  the  imaging  modality  notwithstanding,  the  theoretical  framework 
to  derive  elasticity  using  2-D  and  3-D  motion  models  have  initially  been  developed  by 
members  of  our  group  [24-26].  We  have  also  developed  an  algorithm  to  uniquely 
reconstruct  spatial  maps  of  Young’s  modulus  from  1-dimensional  estimates  of  intern^ 
displacement  and  strain  as  often  done  in  ultrasound  elasticity  imaging.  The  algorithm  has 
been  verified  using  two-dimensional  phantoms  and  a  unique  method  of  ultrasound 
displacement  and  strain  imaging.  Limited  lateral  resolution  of  the  acoustic  imaging 
system  meant  only  one  component  of  the  displacement  vector  was  measured,  but  this  was 
sufficient  given  the  phantom  geometry  to  reconstruct  maps  of  Young’s  modulus. 

In  general,  to  reconstruct  the  tissue-specific  property  of  Young’s  modulus  in 
complex  systems  such  as  the  breast,  the  full  3-D  displacement  vector  must  be  measured 
over  a  3-D  volume.  To  date,  no  group  has  attempted  this  by  any  imaging  modality.  In 
this  project,  we  will  develop  the  techniques  that  can  measure  the  3-D  displacement  vector 
over  any  volume  in  the  object.  These  data  will  be  processed  to  produce  3-D  strain  images 
and  allow  reconstruction  of  Young’s  modulus  over  the  volume.  It  should  be  feasible  to 
translate  these  methods  to  image  the  breast  in  a  clinical  setting. 

METHODS  AND  RESULTS  IN  RELATION  TO  STATEMENT  OF  WORK 

The  goal  of  this  research  program  is  to  develop  a  sensitive  diagnostic  technique 
based  on  quantitative  elasticity  imaging  permitting  surrogate  palpation  of  deep  lying 
breast  lesions.  The  underlying  hypothesis  in  this  project  is  that  measurement  of  elasticity 
in  breast  tissues  provides  unique  information  that  could  increase  detection  and/or 
characterization  of  potentially  malignant  masses  not  accessible  to  manual  palpation.  The 
purpose  of  this  study  is  to  develop  technology  for  acquisition  of  high  resolution  3-D  MRI 
spatial  and  displacement  data  through  the  breast  for  calculation  of  strain  and  elastic 
modulus  images.  To  this  end,  specific  technical  objectives/tasks  were  adopted  and 
summarized  in  the  proposal  Statement  of  Work.  Each  task  is  restated  below  along  with 
descriptions  of  relevant  Methods  and  Progress/Results. 

Technical  Objective  I:  Data  Acquisition  and  Reconstruction 

Task  1:  Month  1-3 

Construction  of  computer  controlled  hydraulic  compression  device  capable  of  producing 
an  incremental  surface  deformation  of  the  mechanical  body.  This  device  will  be 
triggered  by  NMR  imaging  system  and  have  a  simple  timing/displacement  control 
appropriate  for  phantom  studies. 

The  proposed  method  requires  that  the  imaged  object  experience  a  relatively 
minor  externally-applied  deformation  force.  Based  on  signal  optimization  simulations 
(described  below),  differential  deformations  below  10%  should  be  adequate  to  generate 
elasticity  maps  (e.g.,  <lcm  deformation  across  a  10cm  object).  A  drawback,  however,  of 
phase-encode  methods  employed  here  is  that  data  are  acquired  in  many  segments  over  an 
extended  period  of  time  (minutes);  thus,  the  mechanical  deformation  must  be  reproduced 
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for  each  segment.  We  designed  the  acquisition  sequence  to  be  insensitive  to  inevitable, 
ill-defined  motions  (vibrations)  during  deformation  transitions,  but  the  endpoint  of  each 
mechanical  transition  needs  to  be  highly  reproducible.  Clearly,  design  and  fabrication  of 
the  compression  device  is  a  critical  prerequisite  to  phantom  studies.  Development  of  this 
device  dominated  the  first  lOmonths  of  this  project. 

For  electrical  safety  and  handing  ease,  hydraulics  were  abandoned  early  in  the 
design  process  for  pneumatic  drives.  A  high  capacity  >40psi  air  source  piped  into  the 
NMR  laboratory  was  available  to  drive  the  device.  Our  first  automated  deformation 
system  held  the  tissue-mimicking  phantom  between  two  horizontal  plates.  A  pneumatic 
bladder  above  the  top  plate  was  used  to  apply  an  incremental  vertical  deformation  force 
and  mechanical  "stops"  on  guide  posts  were  set  to  define  deformation  limits.  Upon 
release  of  air  pressure  in  the  bladder,  the  phantom  would  passively  recoil  to  its  initial 


Figure  1.  Side-view  schematic  of  pneumatic  deformation  system. 

state.  While  successful  for  initial  proof-of-concept  studies  [27,  28],  we  realized  that  a 
system  based  on  "passive  recoil"  was  inadequate  since  the  recoil  time  itself  was 
dependent  on  viscoelastic  properties  of  the  object  and  the  degree  of  applied  deformation, 
thus  this  timing  would  be  difficult  to  control. 
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Our  second  generation  device  employed  active  deformation  and  active  recoil. 
Machinist  specifications  and  pneumatic  logic  for  the  device  are  shown  in  Figure  1.  Four 
neoprene  bellows  outside  of  the  NMR  coil  (rf-coil)  provide  both  upward  and  downward 
force  to  an  acrylic  plate  on  top  of  the  phantom.  The  top  bellows  are  on  a  common 
pressure  circuit,  likewise  the  bottom  pair  are  pressurized  together  to  yield  uniform 
vertical  displacement  along  the  top  of  phantom.  Solid  acrylic  yolks  at  both  ends  of  the 
device  insure  that  only  the  top  plate  moves  in  response  to  pressurization.  Tests  on  a 
variety  of  bellow  sizes/designs  led  to  current  configuration  of  2.2"  diameter  neoprene 
expansion-joint  boots  with  an  additional  rigid  sleeve  around  each  boot  to  keep  expansion 


Figure  2.  Stimulated-echo  data  acquisition  and  deformation  sequence.  Mechanical 
transitions  occur  during  the  relatively  long  mixing  time  TM.  A  sequence-triggered 
TTL  circuit  controls  pneumatic  relays  to  pressurize  for  active  deformation  and 
recoil. 

forces  vertical.  Further  testing  revealed  that  while  only  one  pair  of  bellows  (top  or 
bottom)  were  pressurized  at  any  one  time,  the  residual  pressure  in  the  opposing  bellows 
greatly  slowed  the  response  of  the  system  to  change  deformation  state.  There  is  a  signal- 
to-noise  advantage  to  rapidly  switch  between  deformation  states.  Quick  exhaust  valves 
for  each  bellow  pair  were  added  to  substantially  improve  system  response  time.  These 
valves  provide  rapid  equalization  to  atmosphere  soon  after  pressure  starts  to  drop  in  the 
bellows. 

Mechanical  stability  tests  of  the  device  were  performed  using  a  slightly  modified 
version  of  the  2D-displacement  encoding  MRI  sequence  illustrated  in  Figure  2. 
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Displacement  encoding  pulses  (Gd)  were  set  to  3.87t/mm  -  15.3;i/mm  sensitivity  along 
horizontal  and  vertical  Erections.  The  deformation  device  was  loaded  with  a  tissue- 
mimicking  phantom,  and  repeated  runs  (128  echoes)  were  acquired  without  phase 
encoding  (i.e.,  G-phase=0)  such  that  phase  variations  in  observed  echoes  would  reflect 
mechanical  instability.  For  example,  assuming  a  displacement  sensitivity  of  3.87i/mm,  a 
measured  phase  instability  of  0.387c  would  suggest  a  mechanical  instability  of  0.1mm. 
The  upshot  of  these  tests  was  that  the  system  is  extremely  stable  and  reproducible  in 
application  of  a  Snun  differential  deformation  (mixing  time,  TM=350msec).  The 
observed  phase  instabilities  were  barely  above  that  allowed  by  NMR  signal  limits.  These 
results  suggest  mechanical  reproducibility  within  <20microns.  Additional  tests  are 
planned  to  determine  stability  at  shorter  mixing  times  of  TM=100  to  250msec. 

Task  2:  Month  1-6 

Development  of  various  phantoms  to  test  NMR  displacement  and  strain  data  acquisition 
and  3-D  elasticity  reconstruction  algorithms. 

Fabrication  of  phantoms  to  model  the  mechanical  and  NMR  properties  of  human 
tissue  is  an  important  step  in  development  and  verification  of  optimized  NMR 
displacement-encoded,  volume  simulated-echo  pulse  sequence,  and  elasticity 
reconstruction  algorithms.  The  phantom  materials,  therefore,  should  closely  resemble 
relevant  tissue  properties.  In  addition,  these  materials  should  be  stable  (i.e.,  material 
properties  should  not  change  over  sufficiently  long  periods  of  time)  and  permit  the  non¬ 
destructive  embedding  of  tissue  specimen  into  the  phantom. 

We  have  previously  developed  tissue-mimicking  phantoms  to  test  NMR  elasticity 
imaging  device.  These  phantoms  were  made  of  two  materials:  a)  polymer  produced  by 
M-F  Manufacturing  Co.,  Inc.  (Fort  Worth,  TX),  and  b)  gelatin  (Sigma-Aldrich  Co.,  3050 
Spruce  Street,  St.  Louis,  MO  63103  USA).  The  first  material,  plastisol,  consists  of  a 
liquid  plastic  combined  with  either  softener  (plasticizer)  or  hardener.  By  varying  the 
proportion  of  these  two  components,  it  is  possible  to  produce  composite  models  of 
desired  elasticity.  The  raw  composite  materids  were  stirred  and  heated  to  approximately 
170°C.  At  that  temperature  the  mixture  was  poured  into  molds  producing  a  tissue- 
mimicking  time-stable  phantom  of  desired  shape  and  elasticity  distribution.  However, 
two  major  deficiencies  remain.  First,  plastisol  does  not  have  desired  tissue  equivalent 
NMR  properties.  Second,  due  to  high  temperature  rise  during  phantom  preparation,  it  is 
impossible  to  produce  tissue-containing  phantoms  using  plastisol. 

Gelatin  phantoms  can  mimic  both  elastic  and  NMR  properties  of  tissue  as  well. 
In  addition,  gelatin  preparation  procedures  allow  fabrication  of  tissue-containing 
phantoms,  but  these  phantoms  cannot  be  kept  over  a  long  period  of  time.  Moreover, 
elastic  properties  of  gelatin  phantoms  are  highly  unstable. 

Several  different  materials  were  further  considered  for  NMR  elasticity  phantoms. 
These  materials  include  Semicosil  921,  Semicosil  905  and  Silgel  612  (Wacker  Silicones 
Co.,  Adrian,  MI  49221),  and  Rhodorsil  RTV  163  (Rhone-Poulenc,  France).  From  all 
tested  materials,  the  Semocosil  921  silicone  gel  appeared  to  mimic  the  tissue  properties 
the  best.  The  elasticity  of  silicone  gel  can  be  simply  controlled  by  mixing  ratio  of  two 
components.  The  composite  material  is  then  poured  into  the  mold  and  cured  at  the  room 
temperature  enabling  fabrication  of  tissue  containing  phantoms.  These  silicone  gels  have 
shown  high  SNR  spin-echo  signal  over  the  range  of  hardness  variations.  In  addition,  the 
mechanical  properties  of  silicone  did  not  change  over  60  days,  as  was  measured  by 
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Instron-type  mechanical  system.  Therefore,  it  is  possible  to  fabricate  tissue-mimicking 
phantoms  of  desired  shape  and  elasticity  distribution  using  Semicosil  921  silicone  gel. 


Task  3:  Month  1-12 

Development  of  the  NMR  displacement-encoded,  volume  simulated-echo  pulse  sequence 
optimized  for  displacement/strain  sensitivity  and  SNR. 

There  are  several  positive  features  in  our  displacement-encoding  stimulated-echo 
sequence  which  are  particularly  well  suited  to  elasticity  imaging.  Difficulties 
encountered  with  “dynamic”  methods  to  document  shear-wave  propagation  are  avoided 
by  the  employed  “static”  displacement  technique.  The  relatively  long  mixing  time  of  the 
stimulated-echo  allows  undesired  dynamic  components  (shear  waves  and  surface 
vibrations)  to  dampen  such  that  static  displacements  are  accurately  encoded. 
Displacement  phase  shifts  are  measured  relative  to  a  reference  file  acquired  with  the 
object  maintained  in  the  same  “State  B”.  This  obviates  the  need  for  image  registration  or 
feature  tracking  algorithms  required  by  other  static  techniques.  In  fact,  local 
displacement  can  be  measured  in  a  featureless  continuum  object.  Moreover,  the  spatial 
resolution  of  displacement  maps  equals  the  resolution  of  the  original  image. 

Despite  these  advantages,  several  parameters  have  yet  to  be  optimized.  In 
particular,  a  good  balance  between  displacement  sensitivity  and  NMR  signal-to-noise 
(SNR)  is  required  since  optimization  of  one  parameter  may  come  at  the  expense  of 
another.  Computer  simulations  were  performed  to  address  these  tradeoffs.  Key 
parameters  that  affect  SNR  and  displacement  sensitivity  include:  displacement  gradient 
amplitude  and  duration;  applied  deformation  amplitude;  image  resolution;  mixing  time, 
TM;  and  tissue  water  diffusion  properties.  Figure  3  illustrates  a  simulation  of  applied 
displacement  gradient  amplitude  and  surface  deformation  (strain)  on  net  SNR.  As  shown, 
a  balance  between  phase  (or,  equivalently,  displacement)  sensitivity  due  to  the  applied 
surface  deformation,  the  magnitude  of  the  encoding  gradient,  and  the  signal  loss  due  to 
diffusion  must  be  obtained  to  maximize  net  SNR.  If  the  combination  of  applied  gradient 
and  surface  deformation  is  too  small,  no  information  about  the  tissue  hardness  will  be 
obtained.  If  the  combination  is  too  large,  intravoxel  dephasing  will  lead  to  a  loss  of  NMR 
signal,  reducing  the  elasticity  image  SNR  as  well  (the  fall-off  shown  with  increasing 
strain  and  gradient).  Finally,  a  displacement  encoding  gradient  that  is  too  large  will  lead 
to  diffusion  losses,  shown  in  Figure  3  as  an  overall  attenuation  as  gradient  is  increased. 
Thus,  these  simulations  suggest  a  net  SNR  gain  via  greater  applied  deformation  with 
reduced  gradient  amplitude.  (Previously,  we  acquired  data  at  4gauss/cm  gradient  strength 
and  <0.02  surface  strain  which  was  well  below  the  optimum  settings).  Mechanical 
reproducibility  in  the  deformation  device  and  patient  tolerance  are  issues  that  need  to  be 
balanced  against  increasing  the  surface  deformation.  The  deformation  system  is  stable  at 
5mm  surface  deformations,  which  we  anticipate  will  be  well  tolerated  by  patients. 

Testing  and  implementation  of  the  3D  displacement-encoding  volume  stimulated- 
echo  sequence  will  begin  once  the  mechanical  deformation  device  is  fully  characterized 
at  shorter  mixing  times.  The  spatial  encoding  segment  of  the  pulse  sequence  (i.e., 
beyond  the  third  rf-pulse)  will  be  extended  in  an  8  -  16  echo  train,  fast-spin-echo  (FSE) 
sequence.  Slice  axis  phase  encoding  will  be  applied  during  the  FSE  echo  train  to  yield 
the  third  spatial  dimension  of  8  -  16  slices.  Addition  of  slice-direction  displacement 
encoding  is  then  straightforward. 


Page  7 


y 


PROPRIETARY  UNPUBLISHED  DATA  Grant  Number:  DAMD 17-97- 1-7079 


SNR  of  Modulus  Reconstruction 


Figure  3.  Relative  signal*to-noise  as  a  function  of  applied  deformation  and 
displacement-encoding  gradient  strength. 

Task  4:  Month  1-12 

Expansion  of  previously  developed  linear  elasticity  reconstruction  algorithms  for 
volumetric  displacement  and  strain  NMR  measurements. 

The  development  and  expansion  of  linear  elasticity  reconstruction  algorithms 
were  primarily  focused  on  the  following  two  aspects:  a)  development  of  boundary 
detection  methods  to  identify  the  regions  of  uniform  elasticity  (Young’s  modulus) 
distribution,  and  b)  development  of  volumetric  linear  elasticity  reconstruction  algorithms 
for  optimized  NMR  displacement-encoded,  volume  simulated-echo  pulse  sequence. 

Boundary  Detection.  The  boundary  detection  is  the  first  step  in  elasticity 
reconstruction,  where  strain  images  are  processed  to  highlight  boundaries  between 
regions  of  different  elastic  modulus.  Following  boundary  detection,  closed  contours  of 
small  elasticity  variations  are  defined.  The  total  dimension  of  the  contour  must  include 
several  resolution  cells  of  the  final  reconstruction.  The  modulus  along  these  contours  can 
be  considered  constant,  thereby  providing  a  boundary  condition  for  complete 
reconstruction  of  the  elasticity  within  this  region  of  interest  (ROI)  based  on  numerical 
solution  of  elasticity  equations. 

The  boundary  detection  algorithm  was  developed  for  specifically  NMR  elasticity 
imaging.  This  algorithm  is  based  on  the  stress  continuity  condition  applicable  for  soft 
tissue  deformations  used  in  these  studies. 
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For  an  inhomogeneous  object,  the  boundary  detection  algorithm  will  give  smoothly 
varying  results  everywhere  except  at  the  boundary  points  where  elasticity  changes.  This 
is  illustrated  in  Figure  4  where  boundary  detection  was  applied  to  the  displacement  and 
strain  data  from  kidney-containing  gel  phantom.  Three  major  regions  of  elastic 
inhomogeneity  can  be  identified  in  Figure  4.  Clearly,  the  gelatin  background  does  not 
exhibit  any  changes  in  Young’s  modulus  throughout  the  entire  imaging  plane.  The  kidney 
parenchyma  appears  as  the  most  elastically  different  material  in  this  phantom.  Finally,  the 
internal  fatty  tissue  in  the  middle  of  the  kidney  is  just  slightly  different  from  the  gelatin. 


Figure  4.  Boundary  detection  algorithm  applied  for  gel-based  phantom 
with  embedded  kidney.  The  longitudinal  view  of  the  kidney  is  depicted  in 
this  image,  where  kidney  parenchyma  and  internal  fatty  regions  are  clearly 
identified. 

In  general,  the  approach  presented  in  Figure  4  produces  a  qualitative  reconstruction,  i.e., 
boundary  detection,  and  will  require  further  more  accurate  Young’s  modulus 
reconstruction  at  every  point  of  the  volume  or  region  of  interest. 

Volumetric  Linear  Elasticity  Reconstruction.  For  a  general,  three-dimensional 
object,  complete  three-dimensional  strain  data  are  required.  However,  to  reconstruct  the 
Young's  modulus  in  a  single  plane,  all  displacement  components  must  only  be  measured 
in  three  closely  positioned  parallel  planes  intersecting  the  volume.  From  measurements 
in  these  three  planes,  coefficients  needed  for  elasticity  reconstruction  in  the  central  plane 
based  on  3-D  equilibrium  equations  can  be  calculated  numerically  with  second  order 
accuracy.  Consequently,  complete  three-dimensional  elasticity  reconstruction  can  be 
done  by  moving  the  imaging  planes  in  parallel.  Note  that  even  a  pair  of  close  parallel 
planes  can  also  be  used  for  numerical  reconstruction  in  a  single  plane,  but  will,  however, 
exhibit  greater  error  since  less  accurate  finite  differences  will  be  used  to  approximate  all 
derivatives. 

Task  5;  Month  10-26 

Development  of  the  nonlinear  elasticity  reconstruction  model  capable  to  process  high- 
strain  NMR  images. 

The  elasticity  reconstruction  problem  can  be  posed  in  number  of  ways  depending 
on  the  experimental  conditions  such  as  amount  of  surface  deformations,  type  of  phantom 
materials,  etc.  Previously,  we  have  developed  a  reconstruction  algorithm  based  on  linear 
theory  of  elasticity,  i.e.,  applicable  for  small  deformations  and  linear  elastic  materials. 
However,  both  mechanical  equilibrium  equations  and  the  displacement-strain  relations 
must  include  high  order  spatial  derivatives  of  the  displacement  if  the  average  deformation 
is  too  large  to  be  described  with  a  linear  model.  An  algorithm  was  developed  to 
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reconstruct  the  elastic  modulus  of  soft  tissue  for  comparatively  large  surface 
deformations  required  for  optimal  SNR  and  sensitivity  strain/displacement  images. 
Numerical  methods  were  developed  to  reduce  error  propagation  in  reconstruction 
algorithms. 

Technical  Objective  II:  Phantom  Studies  on  2T  18cm  Bore  MRI  System 
Task  6:  Month  7-22 

Development  of  gel-  and  rubber-based  phantoms  with  tissue  mimicking  elastic  and  NMR 
relaxation  properties.  Time-stable  phantoms  with  non-palpable  inclusions  of  various 
shapes  and  elasticity  contrast  positioned  at  different  locations  within  the  phantom  will  be 
produced. 

During  the  first  year  of  the  project,  we  have  identified  and  carefully  tested  the 
materials  for  tissue  mimicking  phantoms.  Our  ultimate  goal  is  to  simulate  the  anatomical 
and  geometrical  features  of  the  normal  and  pathologically  transformed  breast  using  these 
materials.  The  models  of  breast  containing  single  or  multiple  inclusions  are  presented  in 
Figure  5.  These  models  were  fabricated  using  plastisol  material  (M-F  Manufacturing 
Co.,  Fort  Worth,  Texas,  USA). 


Figure  5.  Plastisol  based  models  of  breast  with  various  inclusions. 

Our  recent  studies,  however,  suggested  that  silicone  gels  (Semicosil  921,  Wacker 
Silicones  Co.,  Adrian,  Michigan,  USA)  could  better  simulate  the  mechanical  and  NMR 
relaxation  properties  of  the  tissue.  Initial  silicone  gel  based  homogeneous  phantoms  were 
produced  for  mechanical  and  NMR  testing. 

The  following  Technical  Objectives  and  Tasks  will  be  addressed  in  Years  2  and  3  of 
the  project,  but  are  listed  below  for  completeness. 

Task  7:  Month  12-26 

Investigation  of  the  capabilities  of  3-D  NMR  elasticity  imaging,  i.e.,  determine 
sensitivity,  accuracy  and  resolution  of  NMR  displacement  and  strain  images,  and 
reconstructed  elasticity  (Youngis  modulus)  distribution. 

Technical  Objective  III:  Translation  to  1.5T  Human  MRI  System 
Task  8:  Month  20-25 

Design  and  construction  of  the  compression  device  compatible  with  human  breast 
NMR  imaging  system  and  capable  to  produce  wide  range  of  surface  deformations. 
Development  of  sophisticated  circuits  to  control  surface  deformations  and  time 
synchronization  with  human  MRI  system. 
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Task  9:  Month  24-34 

Development  of  appropriate  for  clinical  studies  3-D  displacement  encoded,  volume 
stimulated-echo  pulse  sequence  on  human  MRl  system  with  data  acquisition  time 
within  or  comparable  to  regular  MRI  examination. 

Task  10:  Month  22-34 

Development  of  time  efficient  elasticity  reconstruction  algorithms  more  suitable  for 
clinical  applications. 

Task  11:  Month  20-30 

Estimate  the  influence  of  the  cardiac  and  respiratory  motion  to  elasticity  images  and 
develop  the  approaches  to  reduce  these  artifacts. 

Task  12:  Month  30-36 

Validation  of  clinical  NMR  data  acquisition  and  elasticity  reconstruction  methods  on 
breast  tissue  mimicking  phantoms. 


CONCLUSIONS 

Design,  fabrication,  and  refinement  of  the  deformation  device  essential  for 
phantom  studies  required  more  time  than  originally  planned.  The  current  version  of 
deformation  system,  however,  meets  the  requirements  to  deliver  highly  reproducible 
surface  deformations  of  the  imaged  object.  A  dedicated  TTL  control  circuit  was 
developed  to  adjust  delay  and  duration  of  pulses  that  drive  pneumatic  relays  in  synchrony 
with  in  imaging  sequence.  The  complete  MRI-compatible  hardware,  pneumatic 
components,  and  control  circuitry  are  now  fully  operational.  Pulse  sequence  refinements 
were  investigated  via  simulations  that  indicate  we  were  previously  operating  at  a  sub- 
optimal  regime.  We  have  adjusted  our  acquisition  procedures  accordingly,  and  will  soon 
embark  on  3D  spatial  /  3D  displacement  encoding  phantom  experiments.  Good  progress 
has  been  made  to  refine  3D  elasticity  reconstruction  and  edge  detection  algorithms, 
although  require  additional  phantom  data  is  required  to  further  test  these  strategies.  The 
first  of  these  data  will  be  collected  in  the  next  few  months. 

Significant  effort  was  directed  to  develop  suitable  phantom  materials. 
Requirements  include:  temporally  stable  (over  months-years);  adjustable  mechanical 
properties  to  match  a  range  of  tissues;  moldable  to  simple ^geometries;  tissue-like  NMR 
properties.  The  Semocosil  921  silicone  gel  generally  meets  these  requirements  (with  the 
exception  of  tissue-like  diffusion)  and  will  be  used  in  subsequent  phantom  studies. 
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Elasticity  Reconstructive  Imaging  by  Means  of 
Stimulated  Echo  MRI 

Thomas  L.  Chenevert,  Andrei  R.  Skovoroda,  Matthew  O’Donnell, 
Stanislav  Y.  Emelianov 


A  method  is  introduced  to  measure  internai  mechanicai  dis- 
piacement  and  strain  by  means  of  MRI.  Such  measurements 
are  needed  to  reconstruct  an  image  of  the  elastic  Young’s 
moduius.  A  stimuiated  echo  acquisition  sequence  with  addi- 
tionai  gradient  puises  encodes  internai  dispiacements  in  re¬ 
sponse  to  an  externaiiy  appiied  differentiai  deformation.  The 
sequence  provides  an  accurate  measure  of  static  dispiace- 
ment  by  iimiting  the  mechanicai  transitions  to  the  mixing 
period  of  the  simuiated  echo.  Eiasticity  reconstruction  in- 
voives  definition  of  a  region  of  interest  having  uniform 
Young’s  moduius  along  its  boundary  and  subsequent  soiution 
of  the  discretized  eiasticity  equiiibrium  equations.  Data  acqui¬ 
sition  and  reconstruction  were  performed  on  a  urethane  rub¬ 
ber  phantom  of  known  eiastic  properties  and  an  ex  vivo  ca¬ 
nine  kidney  phantom  using  <2%  differential  deformation. 
Regionai  eiastic  properties  are  weii  represented  on  Young’s 
moduius  images.  The  iong-term  objective  of  this  work  is  to 
provide  a  means  for  remote  paipation  and  elasticity  quantita¬ 
tion  in  deep  tissues  otherwise  inaccessibie  to  manual  palpa¬ 
tion. 

Key  words:  eiastic  Young’s  moduius;  magnetic  resonance 
imaging;  eiastography;  strain  imaging. 

INTRODUCTION 

It  is  well  known  that  tissue  elastic  properties  may  be 
altered  by  tumors.  Young’s  elastic  moduli  may  differ  by 
orders  of  magnitude  in  soft  tissues  in  various  physiologic 
states  (1,  2).  This  finding  is  the  physical  basis  behind 
manual  palpation  used  to  detect  “hard”  masses  (3,  4). 
Indeed,  physical  examination  is  the  first  diagnostic  line 
of  defense  against  breast  cancer,  because  nodule  hard¬ 
ness  raises  suspicion  of  malignancy.  Detection  of  a  new 
breast  mass  by  physical  examination  is  often  sufficient 
for  surgical  excisional  biopsy,  even  when  not  corrobo¬ 
rated  by  other  diagnostic  tests.  Manual  palpation  of  the 
prostate,  superficial  lymph  nodes,  and  abdominal  organs 
are  also  commonly  performed.  Unfortunately,  sensitivity 
of  palpation  is  relatively  poor  within  deep,  dense,  or 
heterogeneous  tissues.  Although  the  touch  of  a  skilled 
interpreter  is  considered  a  powerful  diagnostic  instru¬ 
ment,  most  lesions  detected  by  palpation  tend  to  be  rel¬ 
atively  large  and  superficial. 
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Scientists  are  attempting  to  electronically  extend  the 
touch  of  the  physical  examiner  using  a  variety  of  image- 
based  techniqnes  that  infer  tissue  elasticity.  The  essential 
element  is  measurement  of  internal  motion  and  strain  in 
tissue  structures  experiencing  mechanical  stress.  To 
date,  most  “eiastography”  has  used  ultrasound  to  track 
the  relative  motion  of  targets  by  specular  reflection  (5-7), 
by  Doppler  techniqnes  (8-10),  by  cross-correlation  of 
raw  or  processed  aconstic  echoes  (11, 12),  or  by  tracking 
speckle  patterns  (13-15).  Usually  an  external  static  or 
dynamic  deformation  is  applied  while  internal  displace¬ 
ments  or  propagating  shear  waves  are  documented  by 
imaging. 

MRI  has  also  been  used  to  measure  internal  displace¬ 
ment  and  strain  components  of  the  heart  using  spatial 
magnetization  tagging  (16,  17)  and  phase-based  velocity 
encoding  (18).  Elasticity  reconstruction  of  an  externally 
deformed  phantom  was  demonstrated  using  magnetiza¬ 
tion  tagging,  but  this  method  has  spatial  resolution  lim¬ 
ited  by  the  tagged  grid  size  and  only  measures  2D  motion 
(19).  More  recently,  motion  phase  encoding  by  means  of 
bipolar  gradients  was  used  to  produce  two-dimensional 
(2D)  displacement  and  strain  maps  in  media  mechani¬ 
cally  driven  by  external  forces  (20-23).  Strain  and  dis¬ 
placement  maps  infer  internal  elasticity  but  are  also 
strongly  affected  by  the  applied  deformational  geometry. 
Consequently,  these  maps  do  not  uniquely  reflect  inter¬ 
nal  tissue  properties  (i.e.,  elastic  Young’s  modulus). 
Maps  of  dynamic  strain-wave  propagation,  however,  do 
allow  measurement  of  local  strain  wavelength  or  velocity 
from  which  the  local  elastic  modulus  can  be  derived  (20). 
Shear-wave  attenuation,  interference  from  standing 
waves  off  multiple  reflectors,  and  limited  resolvable 
points  over  the  shear  wavelength  are  potential  drawbacks 
of  this  approach. 

Relative  to  ultrasound,  MRI  has  an  advantage  in  overall 
resolution  and  accuracy  for  multidimensional  displace¬ 
ment  and  strain  measurement  needed  for  elasticity  re¬ 
construction.  Ultrasound  can  accurately  measure  axial 
(i.e.,  along  the  beam  axis)  motion  at  high  spatial  resolu¬ 
tion  (^millimeter),  but  lateral  displacement  is  measured 
at  much  lower  spatial  resolution  defined  by  the  depth- 
dependent  acoustic  heam  width.  The  third  dimension  is 
generally  not  even  considered  given  the  limitations  of 
ultrasound.  Consequently,  reduced  motion  dimensional¬ 
ity  and  overall  low  motion  resolution  of  the  imaging 
system  compromise  the  “elastogram.”  This  shortcoming, 
in  turn,  constrains  the  mechanical  model  used  in  elastic¬ 
ity  reconstruction.  To  date,  only  ID  motion  models  have 
been  applied  to  ultrasound-derived  elasticity  images  of 
tissues  in  vivo  (12,  24).  More  accurate  elasticity  images 
are  achieved  by  properly  controlling  external  deforma- 
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tioiis,  leading  to  2D  elasticity  reconstruction  within  the 
imaging  plane  (25,  26). 

In  this  work  we  present  a  method  to  spatially  encode 
internal  displacement  of  an  object  that  has  undergone  an 
externally  applied  “static”  deformation  with  subsequent 
reconstruction  into  elasticity  maps.  Unlike  dynamic 
techniques  designed  to  estimate  elasticity  from  observa¬ 
tions  related  to  strain-wave  propagation,  static  elasticity 
reconstruction  involves  estimation  of  local  strain  from 
displacement  and  numerical  solution  of  differential  elas¬ 
ticity  equilibrium  equations. 


THEORY  OF  RECONSTRUCTIVE  ELASTICITY 
IMAGING 


The  goal  of  elasticity  imaging  is  to  reconstruct  the  elastic 
modulus  of  a  desired  tissue  region  using  available  mea¬ 
surements  of  displacement  and  strain  components.  In¬ 
deed,  the  mechanical  properties  of  tissue  are  ultimately 
linked  to  the  patterns  of  internal  deformations,  but  the 
deformational  geometry  can  greatly  affect  these  patterns 
as  well.  To  uniquely  image  tissue  elasticity,  the  Young’s 
modulus  must  be  reconstructed  from  estimates  of  inter¬ 
nal  displacement  and  strain. 

In  this  paper,  the  general  approach  to  elasticity  recon¬ 
struction  was  based  on  a  model  of  linear,  elastic,  isotro¬ 
pic,  incompressible  media  (26,  27).  The  key  equations 
and  considerations  are  briefly  presented  here.  A  more 
detailed  description  of  elasticity  reconstruction  is  given 
in  an  earlier  publication  (26). 

In  linear  elasticity,  the  components  of  the  strain  (gy) 
and  stress  (o-ij)  tensors  under  static  deformation  are: 
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1  I  dUi  dUj\ 

2  \ax/  dxj 
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=  -(Ui,j+  Uj.j) 


[1] 


2 

o-ij  =  p8jj  +  -Esi]  [2] 

where  U;  is  a  component  of  the  displacement  vector  U  = 
(Uj,  U2,  U3)  in  Cartesian  coordinates  X  =  (x^,  x^,  X3),  p  is 
the  static  internal  pressure,  Sjj  is  the  Kronecker  delta 
symbol,  and  E  =  E{x^,  X2,  X3)  is  the  Young’s  elastic 
modulus.  Note  in  Eq.  [1],  and  the  entire  paper,  the  lower 
index  after  a  comma  means  differentiation  with  respect 
to  the  corresponding  spatial  coordinate. 

The  static  deformation  of  the  medium  can  be  described 
by  the  equilibrium  condition: 


3 

E  o-i;,;  + /i  =  0,  i=  1,2,3  [3] 

;=i 

where  f,  is  the  body  force  per  unit  volume  acting  in  the  Xj 
direction.  In  addition,  volume  conservation  for  an  in¬ 
compressible  medium  leads  to  the  following  relationship 
between  displacement  and  strain  components 

V  •  [/  =  e-ii  I-  622  +  633  ~  Ui,!  “t  U2,2  +  03,3  =  0  [4] 

Using  Eqs.  [1]  and  [2]  for  stress  and  strain  components, 
and  the  incompressibility  Eq.  [4],  the  equilibrium  condi¬ 


tion  with  eliminated  internal  pressure  p  can  be  rewritten 
in  the  following  form: 

26i2(£'>11  E,22)  +  2(u2>2  n  2e23-^>13 

”  25-12^,23  “  0)12^2)  ^^2 

^12>3-^»3  "I"  3(/2j1  “  fl>2)  ~  0 

2ei3(£’,ii  —  +  2623^)12  2(113,3  “  Wl>l)-^»13 

—  2ei2E,23  +  (V^U3  +  0)13, q)]?,!  ^^>13)2^, 2  [5] 

—  (V^Uq  —  0)q3,3)E,3  +  V^O)q3E  +  3(/3,q  “/qjs)  =  0 

2623(^*22  “  E,23  )  -t-  2613^,32  ^^12^,13  2(U3,3 

—  U2,2)E,23  -I"  Ci>23,iE,i  +  (V^U3  I”  (1)23,2)  E, 2  ~  (V^122 
~  ^ZS’s)  E,3  -t  V^(023E  +  3(/3,2  ~  fzis)  ~  0 


**12  ”  02,1  Ul,2>  **13  ■”  03,1  Ul,3,  **23  ~  ^3>2  ^2<3 

Clearly,  the  elasticity  reconstruction  process  based  on 
Eq.  [5]  requires  accurate  measurements  of  the  displace¬ 
ment  vector,  or,  to  be  more  precise,  requires  accurate 
estimation  of  up  to  third-order  spatial  derivatives  of  the 
displacement.  Equation  [5]  can  also  be  written  in  terms  of 
spatial  derivatives  of  strain  tensor  components.  The 
unique  solution  of  the  system  of  coupled  partial  differ¬ 
ential  equations  [5]  is  determined  by  the  boundary  con¬ 
ditions,  i.e.,  the  elastic  modulus  E(X)  must  be  given  along 
the  boundary  of  the  reconstruction  region  of  interest 
(ROI).  It  should  also  be  noted  here  that  the  analytical 
solution  of  Eq.  [5]  is  not  generally  possible,  and  numer¬ 
ical  methods  must  be  developed  to  solve  this  system  of 
partial  differential  equations. 

Based  on  the  particular  geometry  of  the  phantoms  and 
deformation  system  used  in  these  experiments,  Eq.  [5] 
was  simplified.  A  2D  approximation  of  Eq.  [5]  was  used 
because  the  imaged  plane  was  near  the  center  of  the 


FIG.  1 .  Stimulated-echo  data  acquisition  and  object  deformation 
sequence.  Mechanical  transitions  occur  during  the  iong  mixing 
time  (TM)  such  that  a  static  dispiacement  equilibrium  is  achieved. 
Local  displacement  between  deformation  states  “A”  and  “B”  are 
encoded  by  phase  shift  proportionai  to  “G-displace”  amplitude 
Gc  and  duration  t. 
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phantom.  For  such  a  plane  (arbitrarily  denoted  as  X3  =  0), 
the  displacement  vector  components  do  not  vary  signif¬ 
icantly  as  a  function  of  the  out-of-plane  X3  coordinate, 
and  therefore,  the  “plane  strain  state”  condition  is  appli¬ 
cable.  With  this  condition,  Eq.  [5]  reduces  to  a  single 
nontrivial  equation: 


founded  by  interference  of  the 
primary  shear  wave  with  re¬ 
flected  or  standing  shear 
waves.  To  avoid  this  condi¬ 
tion,  a  “static”  displacement 
encoding  approach  was 
adopted.  It  requires  measure¬ 
ment  of  internal  displacement 
between  two  or  multiple  de¬ 
formations  while  the  object  is 
in  mechanical  eqnilibrium  for 
each  measurement.  A  stimu¬ 
lated  echo  sequence  with  dis¬ 
placement-encoding  gradient 
pulses  is  used  to  achieve  this, 
as  shown  in  Fig.  1.  Mechanical 
transition  from  state  “A”  to 
state  “B”  occurs  during  the 
stimulated  echo  mixing  time, 
TM.  A  relatively  long  mixing 
time  allows  long-lived  elastic 
vibrations  to  dampen  before 
spatial  encoding.  Because  the 
relevant  magnetization  during 
TM  is  longitudinal,  it  is  unaf¬ 
fected  by  potentially  ill-de¬ 
fined  motions  during  the  me¬ 
chanical  transition  period.  As 
a  result,  a  more  accurate  static 
deformation  measurement  is 
achieved.  Also  note  that  pre¬ 
cise  synchronization  of  me¬ 
chanical  and  pulsed  gradient 
events  is  not  critical  as  long  as 
the  mechanical  transition  be¬ 
gins  after  the  second  RF  pulse 
and  is  complete  before  the  third  RF  pulse.  Similarly,  a 
long  delay  in  TE  could  be  used,  but  this  is  done  at  the 
expense  of  signal  lost  to  T2  decay. 

Local  displacement  is  encoded  by  means  of  phase  shift 
governed  by  pulsed-field  gradient  factors. 


FIG.  2.  Urethane  rubber  phantom  containing 
two  hard  inclusions  and  held  between  defor¬ 
mation  plates  on  top  and  bottom  of  the  phan¬ 
tom.  Images  of  stimulated  echo  (a)  magnitude, 
and  phase  shift  proportional  to  (b)  vertical  and 
(c)  horizontal  differential  displacement.  Dis¬ 
placement  sensitivity  of  15.33  w/mm  and  dif¬ 
ferential  displacement  of  =1.4  mm  between 
deformation  plates  was  used. 


{E,ii  —  E,22)(Ui,2  +  U211)  +  2E, 12(112, 2  ~  Ulrl) 

-f  2E, 1(112,11  +  112,22)  ~  2E, 2(111, 11  +  111,22) 

-H  E(U2,iii  -F  U2,221  ~  lll,112  ~  Hi, 222)  +  3(/2>l  “ /I’z)  ~  0 

Under  conditions  in  which  these  assumptions  are  valid, 
only  in-plane  displacement  Ui  and  U2  or  their  derivatives 
are  needed  to  reconstruct  the  modulus  in  the  plane  X3  =  0. 

Static  Displacement  Measurement  by  Means  of 
Stimulated  Echo  MRI 

Shear- wave  propagation  speed  in  soft  tissue  is  1-20  m/s. 
Consequently,  a  shear  wave  imparted  by  a  single-stroke 
or  an  oscillating  deformation  force  may  require  tens  of 
milliseconds  to  traverse  an  object  =100  mm  in  size.  The 
time  for  reflected  waves  to  dampen  may  be  much  longer. 
“Dynamic”  measurements,  which  encode  displacement 
during  shear-wave  propagation,  are  potentially  con¬ 


where  displacement  sensitivity,  <l>d,  is  in  units  of  (radi¬ 
ans/distance).  A  phase  reference  acquisition  is  required 
for  each  displacement  encode  condition  to  remove  pre¬ 
existing  phase  shifts  that  are  unrelated  to  displacement. 
Reference  data  are  acquired  using  the  same  pulse  se¬ 
quence,  including  displacement  encode  gradient  pulses, 
but  with  the  object  maintained  in  state  B.  Note  that  all 
spatial  encoding  occurs  from  the  third  RF  pulse  and 
beyond.  That  is,  the  object  is  in  one  deformation  state 
(state  B)  for  both  displacement  and  phase  reference  ac¬ 
quisitions  during  all  spatial  encoding  segments  of  the 
sequence.  Consequently,  image  registration  or  feature 
tracking  algorithms  (14)  are  not  required  to  estimate  dis¬ 
placement.  Instead,  local  displacement  is  encoded  di¬ 
rectly  by  local  phase  of  a  corrected  dataset,  given  by. 


^cor(l') 


SAf)  Ss(T)* 

lSB(f)| 
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FIG.  3.  Strain  images  calculated  from  first-order  derivatives  of  the  displacement  images  repre¬ 
sented  in  Fig.  2.  Normal  strains  (a)  s,,  and  (b)  S22,  and  (c)  shear  strain  8,2  =  £21  reflect  internal 
elastic  properties  and  the  externally  applied  deformation  field.  The  plane  strain  state  assumption 
and  phantom  incompressibility  suggest  ^  -822  as  is  supported  by  the  relatively  featureless 
map  of  81.,  +  822  in  (d). 


where  and  Sg  are  the  data  acquired  with  the  object 
initially  in  state  A  and  state  B,  respectively.  The  un¬ 
wrapped  phase  of  the  corrected  dataset  and  Eq.  [7]  pro¬ 
vides  a  local  measure  of  displacement,  U,  by  means  of; 

<p(f)  =  4>d  •  [r^  -  ffl]  =  •  U{f)  [9] 

Most  sources  of  phase  error,  such  as  static  field  inho¬ 
mogeneity,  tend  to  be  slowly  varying  functions  of  posi¬ 
tion.  Therefore,  the  phase  reference  datasets  may  be  ac¬ 
quired  at  relatively  low  spatial  resolution  to  reduce  scan 
time. 

METHODS 

Data  Acquisition 

Elasticity  imaging  was  performed  on  two  phantoms.  One 
was  an  85-mm  diameter  cylindrical  urethane  rubber 
phantom  containing  two  8-mm  cylinders  of  hard  mate¬ 
rial.  Previously,  the  ratio  of  Young’s  modulus  between 
the  inclusion  and  background  material  was  measured 


10.5  ±  1.5  for  30%  surface  de¬ 
formation  (19).  In  the  present 
study,  smaller  surface  defor¬ 
mation  was  applied,  and  con¬ 
sequently,  the  elasticity  con¬ 
trast  within  this  phantom 
should  be  less  than  the  previ¬ 
ously  measured  ratio.  Whereas 
the  mechanical  properties  of 
the  rubber  phantom  mimic 
soft  tissue,  the  phantom  had 
inherently  poor  NMR  signal.  A 
more  tissue-equivalent  phan¬ 
tom  in  terms  of  NMR  and  me¬ 
chanical  properties  was 
achieved  by  embedding  a  fresh 
canine  dog  kidney  (<24  h  ex 
vivo]  into  a  130  mm  X  105 
mm  X  75  mm  block  of  5%  gel¬ 
atin.  One  hour  before  MRI,  10 
ml  of  5%  glutaraldehyde  solu¬ 
tion  was  injected  into  kidney 
parenchyma  to  create  a  hard 
lesion(s). 

Phantoms  were  held  se¬ 
curely  in  place  under  moder¬ 
ate  preload  pressure  by  two 
parallel  acrylic  plates.  When 
pressure  to  the  top  plate  was 
released,  the  phantom  recoiled 
vertically  an  amount  con¬ 
strained  by  physical  stops. 
Maximum  vertical  displace¬ 
ment  was  <1.5  mm,  which 
represented  <2%  differential 
between  state  “A”  (greater  de¬ 
formation)  and  state  “B”  (less 
deformation).  Deformation 
was  actuated  pneumatically 
by  an  air-filled  bladder  on  top 
of  the  phantom  holder.  Pneu¬ 
matic  pressure  was  stepped  by  a  remote  solenoid  valve 
with  timing  controlled  by  an  external  transistor  transistor 
logic  (TTL)  gate  circuit  triggered  by  the  pulse  sequence. 

Displacement  encoding  gradient  pulse  duration,  t  = 
4.5  msec,  and  amplitude,  Gd  =  40  mT/m,  provided  a 
displacement  sensitivity  of  3)^  =  15.33  Tr/mm  by  means 
of  Eq.  [7].  The  displacement  encoding  direction  was  al¬ 
ternated  each  pulse  repetition  between  vertical  and  hor¬ 
izontal.  For  the  urethane  rubber  phantom,  acquisition 
pareuneters  were  TR  =  1.3  s,  TM  =  350  msec,  TE  =  50 
msec,  128  X  128  matrix,  four  signal  averages,  100  mm 
field  of  view,  and  10-mm  section  thickness.  An  addi¬ 
tional  128  X  32  matrix  acquisition  was  collected  for 
phase-reference  correction  of  vertical  and  horizontal  en¬ 
coded  data.  The  kidney  phantom  acquisition  parameters 
were  TR  =  1  s,  TM  =  200  msec,  TE  =  76  msec,  256  X  256 
matrix,  two  signal  averages,  150  mm  field  of  view,  and 
5-mm  section  thickness;  with  a  256  X  32  dataset  acquired 
for  phase  correction.  All  experiments  were  performed  on 
a  2  T  18-cm  bore  MRI  system  (Bruker,  formerly  GE  NMR 
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FIG.  4.  (a)  Reconstructed  map  of  Young’s  modulus  for  the  urethane 
rubber  phantom  within  a  63  x  50  mm  region.  The  boundary  of  the 
region  is  defined  to  have  a  Young’s  moduius  =  1 .  Reiative  Young’s 
moduii  for  inciusions  and  background  materiai  are  shown  in  (b). 

Instruments),  using  a  150-mm  transmit/receive  birdcage 
coil. 

Data  Processing 

Time-domain  data  were  transferred  for  off-line  process¬ 
ing  as  follows.  Phase  reference  datasets  were  zero-filled 
and  2D  Fourier  transformed  to  a  128  X  128  or  256  X  256 
matrix  for  phase  correction  by  means  of  Eq.  [7].  The 
resulting  phase  maps  were  used  to  estimate  the  spatial 
derivatives  of  the  in-plane  displacements  necessary  for 
elasticity  reconstruction  (6).  Note  phase  unwrapping  is 
not  strictly  required  because  only  phase  derivatives  are 
used.  Assuming  the  displacement  fields  are  continuous, 
resulting  in  small  differential  displacement  at  any  pixel 
compared  with  the  total  displacement,  the  differential 
displacement  between  two  neighboring  pixels  was  di¬ 
rectly  computed  from  the  angle  of  the  complex  multipli¬ 
cation  of  each  pixel  with  the  conjugate  of  the  neighboring 
pixel,  then  scaled  by  l/<hd' 

Solving  Eq.  [6]  for  unknown  E  (xl,  x2)  performed  the 
elasticity  reconstruction,  i.e.,  reconstruction  of  the  spa¬ 
tial  distribution  of  elastic  Young’s  modulus.  As  was 
noted  previously,  the  unique  solution  of  a  boundary 
value  problem  (Eq.  [5]  or  [6])  is  determined  by  the  bound¬ 


ary  conditions.  Therefore,  a  rectangular  ROI  was  identi¬ 
fied  within  the  imaging  planes  for  both  phantoms.  For 
the  phantom  with  two  hard  inclusions,  the  ROI  was  a 
region  of  63  X  50  mm  positioned  approximately  in  the 
center  of  the  phantom  and  included  both  inclusions.  For 
the  canine  kidney  phantom,  the  rectangular  94  X  51  mm 
ROI  included  the  whole  kidney  cross-section.  In  both 
cases,  the  Young’s  modulus  value  along  the  ROI  boundary 
was  set  to  “one”  resulting  in  reconstruction  of  relative 
Young’s  modulus.  More  detailed  analysis  and  discussion  of 
defining  the  ROI  was  considered  previously  (26). 

Elasticity  reconstruction  Eqs.  [5]  and  [6]  assume  that 
spatial  derivatives  of  the  Young’s  modulus  are  continu¬ 
ous  functions.  To  ensure  continuous  elasticity  distribu¬ 
tion,  the  spatial  derivatives  of  the  displacement  were 
low-pass  filtered  before  elasticity  reconstruction,  result¬ 
ing  in  mild  spatial  resolution  reduction. 

After  defining  the  boundary  conditions,  Eq.  [6]  was 
discretized  over  the  ROI  with  the  same  grid  spacing  as 
the  MR  images,  where  all  spatial  derivatives  of  the  dis- 
placement/strain  (i.e.,  coefficients  in  Eq.  [6]  for  unknown 
Young’s  modulus  distribution)  were  approximated , by 
finite  differences.  The  linear  set  of  equations  resulting 
from  discretization  of  Eq.  [6]  was  solved  iteratively, 
where  the  error  in  each  step  was  estimated  by  averaging 
the  left-hand  side  of  Eq.  [6]  over  the  ROI  using  the  current 
estimate  of  the  elasticity  distribution.  From  step  to  step, 
the  Young’s  modulus  distribution  was  updated  based  on 
the  changes  in  the  average  error. 

RESULTS 

Magnitude  and  corrected  phase  images  of  the  urethane 
rubber  phantom  are  shown  in  Fig.  2.  Given  that  <E>d  = 
15.33  ir/mm,  the  number  of  2tt  phase  bands  in  Fig.  2b 
indicates  that  the  vertical  excursion  of  the  phantom  top 
relative  to  the  bottom  was  »=1.4  mm;  similarly,  the  rela¬ 
tive  lateral  displacement  between  left  and  right  edges  of 
the  phantom  was  '“0.8  mm  (Fig.  2c).  Reduced  phase 
slopes  in  the  regions  of  the  hard  inclusions  are  clearly 
visible  on  the  phase  images.  Normal  strain,  and  822, 
and  shear  strain,  —  E21,  maps  are  illustrated  in  Figs. 
3a-3c,  respectively.  The  observed  contrast  reversal  be¬ 
tween  Sij  and  E22  is  a  result  of  phantom  incompressibil¬ 
ity  (like  soft  tissue),  which  yields  =  — e22  assuming 
negligible  out-of-plane  strain.  Consequently,  the  sum  (e^ 
+  S22)  is  relatively  “flat”  as  shown  at  equivalent  gray¬ 
scale  settings  in  Fig.  3d.  Also  note,  whereas  the  strain 
maps  clearly  exhibit  object-specific  detail  (i.e.,  inclu¬ 
sions),  features  related  to  the  applied  external  deforma¬ 
tion  are  quite  conspicuous.  This  fact  demonstrates  why 
an  elasticity  reconstruction  is  needed.  The  Young’s  mod¬ 
ulus  was  reconstructed  for  a  63  X  50  mm  region  as 
presented  in  Fig.  4a.  The  boundary  of  the  elasticity  re¬ 
construction  area  was  defined  to  have  value  “1.”  The 
relative  elastic  moduli  for  select  regions  are  indicated  in 
Fig.  4b  and  are  consistent  with  the  known  elasticity  of 
these  materials  (19). 

The  canine  kidney  phantom  is  shown  in  Fig.  5  as 
magnitude  (Fig.  5a),  vertical  (Fig.  5b),  and  horizontal 
(Fig.  5c)  phase  shift  images.  It  is  apparent  from  the  ver¬ 
tical  phase  image  (Fig.  5b)  that  although  '«1.24-mm  rel- 
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FIG.  5.  Stimulated  echo  images  of  canine  kid¬ 
ney  imbedded  in  gelatin,  (a)  Magnitude  image, 
and  (b)  vertical  and  (c)  horizontal  displacement 
phase  images  in  which  <1.3  mm  of  verticai 
differential  deformation  and  15.33  ir/mm  dis¬ 
placement  sensitivity  were  applied. 


C 

ative  displacement  spans  the  full  phantom,  only  *=0.26 
mm  of  it  is  within  the  kidney.  Consequently,  there  is  a 
high  concentration  of  strain  near  the  gel-kidney  interface, 
as  is  clear  on  the  strain  maps  (Fig.  6).  As  noted  before, 
strain  images  exhibit  contrast  related  to  a  combination  of 
internal  structure  and  the  externally  applied  deforma¬ 
tion.  Elasticity  reconstruction,  however,  reduces  the  am¬ 
biguity  and  exhibits  contrast  dominated  by  internal  elas¬ 
tic  properties  (Fig.  7a).  It  is  also  encouraging  to  note  that 
whereas  there  was  only  moderate  strain  contrast  within 
the  kidney,  elasticity  contrast  within  renal  parenchyma 
and  central  sinus  are  well  distinguished  on  the  Young’s 
modulus  image.  Moreover,  the  site  of  glutaraldehyde  in¬ 
jection  (top-right  quadrant  of  images)  exhibited  the  high¬ 
est  relative  Young’s  modulus.  Approximately  20  h  after 
MRI,  the  kidney  phantom  was  sliced  at  a  plane  corre¬ 
sponding  to  that  studied  by  MRI.  An  optical  image  of  this 
slice  is  shown  in  Fig.  7c.  The  freshly  cut  smface  was  pal¬ 
pated  such  that  areas  of  relatively  “hard”  parenchyma 
could  be  noted.  Arrows  in  Fig.  7c  mark  the  most  conspic¬ 
uous  areas  of  hardness;  the  largest  area  corresponds  to  the 
high  Young’s  modulus  region  in  the  upper-right  quadrant 
of  the  kidney. 

DISCUSSION 

In  this  work  we  introduce  a  method  to  image  and  quan¬ 
tify  internal  elastic  properties  of  an  object  by  means  of 
displacement-sensitive  MRI  with  associated  elasticity  re¬ 
construction.  The  data  acquisition  segment  employs  gra¬ 
dient  pulses  to  encode  internal  displacement  by  means  of 


phase  using  a  stimulated  echo 
sequence.  Internal  displace¬ 
ments  occur  in  response  to  an 
external  deformation  force 
synchronized  to  the  acquisi¬ 
tion  sequence.  By  timing  con¬ 
trol,  mechanical  motion  oc¬ 
curs  while  the  relevant 
magnetization  is  longitudinal. 
The  stimulated  echo  allows 
extension  of  the  mechanical 
transition  period  to  avoid  po¬ 
tentially  long-lived  or  ill-de¬ 
fined  oscillations  within  the 
object  such  that  an  estimate  of 
“static”  displacement  is 
achieved.  Image  registration  or 
feature  tracking  (14)  is  not  re¬ 
quired  because  the  object  is  in 
one  deformation  state  for  all 
spatial  encoding,  fy  relaxation 
and  diffusion,  which  erode 
signal,  ultimately  set  practical 
limits  on  this  period.  In  these 
experiments,  200  to  350  msec 
was  sufficient  to  allow  static 
displacement  measurement  of 
rubber  and  gelatin/tissue 
phantoms  using  a  simple  air- 
bladder  pneumatic  system. 
This  arrangement  required  the 
object  to  “passively”  recoil  during  the  TM  period. 
Clearly,  a  faster  deformation  system  can  be  built  using 
external  forces  to  “actively”  deform  the  object  during  the 
stimulated  echo  mixing  period.  A  shorter  mixing  time 
would  then  be  used  yielding  a  larger  signal. 

To  date,  two  approaches  are  present  in  elasticity  im¬ 
aging:  static  reconstructive  elasticity  imaging  (25,  26)  and 
dynamic  shear-wave  elasticity  imaging  (20,  21,  23).  In 
both,  an  external  static  or  dynamic  deformation  is  ap¬ 
plied  while  the  resulting  displacement/strain  or  propa¬ 
gating  shear  wave  is  detected  using  an  imaging  modality. 
In  reconstructive  elasticity  imaging,  the  elasticity  distri¬ 
bution  must  be  reconstructed  from  static  displacement 
and  strain  images.  The  ability  to  control  the  internal 
deformation  pattern  by  varying  the  externally  applied 
load,  and  high  SNR  displacement  and  strain  estimates 
are  the  benefits  of  this  method,  although  numerical  re¬ 
construction  algorithms  are  required.  Wherein  the  mod¬ 
eled  assumptions  are  valid,  these  algorithms  exist.  For 
more  general  applications,  they  must  be  refined  further. 
In  shear-wave  elasticity  imaging,  local  shear  wavelength 
measurements  allow  direct  and  simple  calculation  of  the 
shear  elastic  modulus.  However,  the  interference  of  shear 
waves  reflected  from  any  elasticity  inhomogeneities 
within  the  tissue,  along  with  attenuation  of  shear  waves, 
and  conversion  between  shear  and  bulk  waves  are  chal¬ 
lenges  of  this  method. 

In  these  experiments  a  displacement  sensitivity  of  <l>d 
=  15.33  Tr/mm  was  achieved  using  moderate  gradient 
factors.  The  ultimate  quality  of  Young’s  modulus  recon¬ 
struction  depends  on  the  induced  phase  shift,  equal  to 
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FIG.  6.  Strain  maps  of  the  kidney  phantom. 
Normal  strains  (a)  e,,  and  (b)  and  (c)  shear 
strain  6,2  illustrate  fairly  high  concentration  of 
strain  near  the  gel-kidney  interface. 
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the  product  of  and  local  displacement.  There  is  a 
limit,  however,  to  the  advantages  gained  by  increasing 
phase  shift.  As  spatial  phase  gradients  become  large,  the 
phase  distribution  within  a  given  voxel  reduces  signal 
amplitude.  Assuming  a  linear  phase  distribution  of  range 
j3  with  a  voxel,  the  signal  modulation  function  is  given  by 
sinc(j8).  For  illustration,  consider  the  vertical  phase  ex¬ 
cursion  of  21  17  observed  across  the  80-mm  phantom  in 
Fig.  2b.  If  the  vertical  phase  excursion  was  evenly  dis¬ 
tributed  across  voxels,  the  phase  range  within  each 
0.78-mm  voxel  would  be  approximately  0.2  ir.  This  im¬ 
plies  a  signal  reduction  factor  of  sinc(0.2  it)  =  0.94  (i.e., 
signal  loss  of  6%).  Clearly  phase  gradients  can  be  more 
concentrated  depending  on  object  geometry,  elastic  het¬ 
erogeneity,  and  deformation  geometry.  This  concentra¬ 
tion  can  lead  to  regions  of  significant  signal  loss  in  a 
manner  analogous  to  flow  dephasing  in  conventional 
MRI.  Under  such  conditions,  a  smaller  voxel  size  can 
(paradoxically)  yield  higher  signal.  In  addition,  inspec¬ 
tion  for  significant  signal  loss  within  the  displacement- 
sensitive  magnitude  image  can  identify  high-strain  re¬ 
gions  near  soft/hard  interfaces  of  a  lesion. 

Water  diffusion  in  the  presence  of  displacement  en¬ 
coding  gradients  is  another  source  of  signal  attenuation. 
We  estimate  the  signal  reduction  factor  for  freely  diffus¬ 
ing  water  was  ^0.24  in  these  experiments  (i.e.,  76% 
signal  lost  to  diffusion  effects).  Diffusion  effects  are  less¬ 
ened  by  reducing  4>jj,  or  alternatively  by  shortening  TM 
without  affecting  In  either  case,  diffusion  effects  are 
assumed  independent  of  the  deformation  state,  and 


therefore  are  ignored  in  the 
elasticity  reconstruction.  Be¬ 
cause  displacement  phase 
shift  is  the  product  of  local 
displacement  and  the  se¬ 
lection  of  is  somewhat  ar¬ 
bitrary  as  long  as  the  applied 
differential  deformation  is  ad¬ 
equate  for  elasticity  recon¬ 
struction.  In  these  preliminary 
experiments,  the  differential 
deformation  was  <1.5  mm 
across  the  imaged  object.  For 
multi-step  acquisitions,  as 
done  here,  good  reproducibil¬ 
ity  of  deformation  is  essential. 
Significant  variation  in  defor¬ 
mation  magnitude  over  the  ac¬ 
quisition  will  lead  to  phase  in¬ 
stability,  motion-like  artifact 
in  base  images,  and  errors  that 
propagate  through  the  elastic¬ 
ity  reconstruction.  It  is  a  minor 
technical  challenge  to  achieve 
relatively  high  displacement 
reproducibility  in  the  defor¬ 
mation  apparatus,  hreproduc- 
ible  motions  that  originate 
within  the  imaged  object, 
however,  can  be  problematic 
and  are  analogous  to  undes¬ 
ired  physiologic  motion  arti¬ 
facts  in  in  vivo  diffusion  MRI.  Fortunately,  unlike  diffu¬ 
sion  and  physiologic  motions,  the  targeted  motion  in 
elasticity  imaging  is  externally  driven.  As  such,  the  dis¬ 
placement  amplitude  in  response  to  an  external  differen¬ 
tial  deformation  can  be  significantly  greater  than  irrepro- 
ducible  or  asynchronous  displacement.  For  many  in  vivo 
applications  including  the  breast,  increasing  the  differential 
deformation  severalfold  relative  to  that  applied  in  these 
phantom  studies  can  reduce  motion  artifact.  Gradient  fac¬ 
tors  and  would  be  reduced  accordingly,  which  would 
yield  the  added  benefit  of  increased  signal  otherwise  lost  to 
diffusion  effects. 

In  practice,  the  definition  of  a  closed  contom  of  constant 
Yoimg’s  modulus  within  the  tissue  can  be  a  challenge.  A 
hybrid  procedure  can  be  used  as  detailed  elsewhere  (26) 
and  summarized  as  follows.  The  strain  images  are  first 
processed  to  highlight  bmmdaries  between  regions  of  dif¬ 
ferent  elastic  modulus.  This  procedure  is  based  on  the 
stress  continuity  property  of  continuous  media  such  as 
tissue  and  can  define  regions  of  very  small  modulus  varia¬ 
tions.  After  this  boundary  detection,  closed  contours  of 
small  elasticity  variations  are  defined.  The  modulus  along 
the  contours  is  considered  constant,  thereby  providing  the 
boundary  condition  for  complete  reconstruction  of  the  elas¬ 
ticity  within  the  region  of  interest  based  on  numerical  so¬ 
lution  of  Eq.  [6].  The  elasticity  distribution  reconstructed  in 
this  way  is  the  modulus  relative  to  the  modulus  along  the 
boundary.  For  breast  elasticity  imaging,  it  is  anticipated 
that  such  a  contom  can  be  defined  a  priori  within  the 
subcutaneous  fat  that  surrounds  breast  parenchyma.  As- 
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FIG.  7.  (a)  Reconstructed  Young’s  modulus  Image  within  a  rectangular  94  x  51  mm  region 
encompassing  the  kidney,  (b)  Relative  Young's  modulus  for  the  scribed  ROIs  indicate  high  elastic 
modulus  at  the  glutaraldehyde  injection  site  in  the  upper-right  quadrant  of  the  kidney  parenchyma, 
(c)  The  optical  image  of  the  kidney  phantom  approximateiy  20  h  after  MRI.  Areas  marked  by  arrows 
were  noticeabiy  harder  as  assessed  by  sense  of  manuai  touch. 


range  to  include  deep  lying  le¬ 
sions.  One  application  would 
be  measurement  of  elasticity 
in  breast  tissue  not  accessible 
to  manual  palpation.  In  situ 
studies  of  Young’s  elastic 
modulus  performed  on  sam¬ 
ples  of  breast  tissue  indicate 
that  there  is  a  large  difference 
in  elastic  modulus  between 
normal  and  pathologically 
transformed  breast  tissues. 
Others  have  analyzed  the 
Young’s  modulus  differences 
between  different  soft  tissues 
and  have  found  1-2  orders  of 
magnitude  difference  in 
Young’s  elastic  moduli  of  a  tis¬ 
sue  in  different  physiologic 
states  (1).  If  elastic  changes 
predate  formation  of  calcifica¬ 
tions,  elasticity  imaging  could 
potentially  increase  detection 
and/or  characterization  of  ma¬ 
lignant  breast  masses  and  thus 
be  an  important  addition  to  ex¬ 
isting  clinical  diagnostic  tools. 
Practical  issues  such  as  the  rel¬ 
atively  high  cost  of  MRI  may 
hinder  use  of  this  approach  as  a 
screening  test.  Nevertheless,  ad¬ 
ditional  work  to  define  the  role 
of  this  technique  as  a  primary 
diagnostic  tool  or  supplemental 
problem-solving  modality  in 
the  management  of  soft-tissue 
disease  is  well  justified. 


suming  the  Young’s  modulus  of  the  fat  boundary  is  con¬ 
stant,  the  relative  Young’s  modulus  image  of  parenchyma 
can  be  reconstructed.  Alternatively,  the  breast  can  be  sirr- 
rormded  by  a  high-signal  cuff  of  known  elastic  modulus 
and  imaged.  An  image  of  absolute  Young’s  modulus  can 
then  be  reconstructed  if  the  boundary  contour  is  defined 
within  the  cuff  material. 

Some  artifacts  present  in  elasticity  images  (Figs.  4  and  7) 
are  due  to  violation  of  the  plane  strain  state  approximation 
in  these  experiments.  Indeed,  if  a  plane  strain  state  is  not 
present,  the  reconstruction  based  on  Eq.  [6]  will  be  in  error. 
The  elasticity  reconstruction,  however,  does  not  have  to  be 
limited  by  a  plane  strain  assumption  if  all  3D  components 
of  the  displacement  vector  are  available.  Fortrmately,  such 
information  would  be  available  using  3D  displacement  en¬ 
coding  within  a  volumetric  imaging  sequence.  The  issue  of 
long  scan  time  could  be  resolved  by  incorporating  echo- 
planar  imaging  or  fast-spin-echo  segments  for  spatial  en¬ 
coding.  Correspondingly,  an  elasticity  reconstruction  based 
on  Eq.  [5]  would  be  applied  to  produce  volumetric  elastic¬ 
ity  maps. 

The  long-range  goal  of  quantitative  elasticity  imaging 
is  to  provide  remote  palpation  thus  expanding  its  limited 
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